High field magnetospectroscopy to probe the 1.4 eV Ni color center in diamond 
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A magneto-optical study of the 1.4 eV Ni color center in boron- free synthetic diamond, grown at 
high pressure and high temperature, has been performed in magnetic fields up to 56 T. The data is 
interpreted using the effective spin Hamiltonian of Nazare, Nevers and Davies [Phys. Rev. B 43, 
14196 (1991)] for interstitial Ni^ with the electronic configuration Sd^ and effective spin S — 1/2. 
Our results unequivocally demonstrate the trigonal symmetry of the defect which preferentially 
aligns along the [111] growth direction on the (111) face, but reveal the shortcomings of the crystal 
field model for this particular defect. 

PACS numbers: 71.55.-i, 78.55.-m 



I. INTRODUCTION 

Diamond, as a material, has attracted a lot of atten- 
tion due to its unique physical properties; it is the hard- 
est known material with high thermal conductivity and 
a 5.5 eV wide electronic band gap. The large Debye 
temperature of diamond reduces the interaction between 
impurities and the lattice leading to almost atomic like 
optical emission spectra of defects with extremely nar- 
row lines. The so-called color centers, various transition 
metal- nitrogen/ vacancy complexes in diamond, can act 
as single photon sources, capable of photostable opera- 
tion at room temperature,^^ with possible applications 
in quantum information processing^. Moreover, the ni- 
trogen/vacancy (NV) center in diamond has been used to 
image a single electronic spin using nanoscale magnetom- 
etry.^'^ The NV center has also been used to produce di- 
amond based light emitting diodes.^ These developments 
have led to a renewed interest in the optical properties 
of color centers in diamond. 

Macroscopic synthetic diamond crystals are mainly 
grown by chemical-vapor deposition (CVD) methods or 
by the high-pressure high-temperature (HPHT) method, 
the latter giving mm^ size bulk crystals. In both growth 
techniques there are only a few impurities which can enter 
into the diamond structure. Using the HPHT method, 
the incorporation of cobalt (Co)^ or nickel (Ni)^^'^^ has 
been achieved in significant amounts. Incorporating tran- 
sition metals into diamond is of interest for applica- 
tions such as single photon emitters^"^ or spintronics 
in analogy to diluted magnetic semiconductors (DMS). 
The Curie temperature of DMS has been predicted to 
scale with the inverse cube of the lattice constant of the 
host matrix^^. Diamond has the smallest lattice constant 
(a = 3. 56 A) of all semiconductors, making it an excellent 



candidate for ferromagnetic ordering above room temper- 
ature. For all these reasons, a thorough understanding 
of transition metal complexes in diamond has become es- 
sential to further develop diamond related technologies. 

Despite the numerous potential applications of the 
nickel color center in diamond, its exact crystallographic 
site and electronic properties are still under debate. The 
NIRIM-2 electron spin resonance (ESR) line,^^ has been 
identified with the 1.4 eV doublet of zero phonon lines 
(ZPLs) seen in optical studies, ^^^^ and unambiguously 
attributed to a nickel containing center. ^^'^^ The ESR 
results suggested that the nickel is incorporated intersti- 
tially in a single positively charge state 3(i^, effective spin 
6* = 1/2 with trigonal symmetry and a strong trigonal 
distortion due to the presence of an additional impurity 
or vacancy near by^^. The trigonal symmetry was con- 
firmed by optical studies^^ under uniaxial stress and mag- 
netic fields up to 6 T. The agreement between the ESR 
and magneto-optical data has lead some authors to pro- 
pose that the Zeeman splitting of the 1.4 eV line be used 
as a pulsed magnetic field calibration probe. ^^ However, 
recent studies suggested an alternative complex of nickel 
with boron, ^^ or even isolated interstitial Ni.^^ Recent 
theoretical work also contradicts the hypothesis of either 
isolated interstitial Ni or interstitial Ni complex with ei- 
ther a vacancy or an impurity; the calculations predict 
that interstitial Ni with trigonal symmetry is unstable. ^^ 

In order to further elucidate the nature of this de- 
fect, we have performed a magneto-optical study of the 
Ni color center in synthetic diamond which is charac- 
terized by two zero-phonon lines (ZPLs) at 1.401 and 
1.404 eV. Measurements in high magnetic fields unequiv- 
ocally demonstrate the trigonal symmetry of the defect 
which is mostly incorporated on the (111) faces with its 
trigonal axis preferentially aligned along the growth di- 



rection. For incorporation on the (001) face there is no 
preferential ahgnment of the defect axis along one of the 
four < 111 > directions. These results are extremely well 
described by the effective spin Hamiltonian for a trigonal 
defect proposed in the seminal paper of Nazare, Nev- 
ers and Davies^^ (hereafter referred to as NND) using a 
single parameter set for both the (111) and (001) faces. 
However, in certain regions of the (111) face the trigo- 
nal axis of the defect is not preferentially aligned. Under 
such conditions, the correct parameters of the effective 
spin Hamiltonian are significantly different. 



II. SAMPLE CHARACTERIZATION 

Two diamond crystals containing nickel and nitrogen 
have been studied. The crystals were grown in nickel 
solvent in similar high pressure and high temperature 
conditions, except that for one of the crystals a nitrogen 
getter (Ti) was added to the solvent in order to reduce 
the nitrogen incorporation. Crystal KA1970 (Sample A), 
grown with the nitrogen getter has a light green color, 
and, KA1153 (Sample B) grown without the nitrogen 
getter has a yellow/brown color. As there was no specific 
boron contamination of the growth environment, the con- 
centration of this impurity is expected to be negligible in 
the grown crystals. Both single crystals have mainly wide 
(111) growth sectors, with smaller (001) growth sectors 
terminated by square faces. The type lb diamond seed 
crystals were not removed. According to the growth con- 
ditions, the nickel concentration is expected to be around 
10^^ cm~^ in both crystals, and the nitrogen concentra- 
tion around 10^^ cm~^ in Sample A and 10^^ cm~^ in 
Sample B. Due to the dependence of the impurity incor- 
poration rate on the crystallographic orientation during 
growth, the impurity concentration is expected to de- 
pend on the growth sector. Crystals were characterized 
by cathodoluminescence at T = 5 K with a 30 kV e-beam 
and magnetometry using a SQUID magnetometer in the 
temperature range 2 to 300 K and in a to 5 T magnetic 
field. 



A. Cathodoluminescence 

For Sample B, grown without the nitrogen getter, 
the cathodoluminescence spectra recorded on (100) and 
(111) growth sectors exhibits mainly the 1.40 eV center 
(Fig. 1), which can be attributed unambiguously to a de- 
fect containing a single nickel atom^^. Nitrogen related 
centers H3 (2.46 eV, attributed to nitrogen VNV com- 
plex) and 2.156 eV (attributed to neutral NV complex) 
also appeared, but with stronger H3 signal on the (100) 
sectors and stronger 2.156 eV signal on (111) sectors. A 
nickel related peak at 2.56 eV was also observed on the 
(111) growth sector, suggesting a stronger incorporation 
of Ni in the (111) growth sectors, in agreement with pre- 
vious studies. Note that no signal corresponding to W8 
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FIG. 1. (Color on-line) Cathodoluminescence spectra of Sam- 
ple A (grown with N getter) and Sample B (grown without 
N getter) recorded on a (100) and (111) growth sectors under 
similar conditions. 



centers, related to substitutional Ni~ were observed. ^^'^^ 
For Sample A, grown with the nitrogen getter, the inten- 
sity of the 1.40 eV lines is significantly larger compared 
to Sample B. Spectra recorded on the (111) growth sector 
reveal the 2.56 eV Ni related center, with very weak N 
related signals. This suggests a much stronger incorpo- 
ration of Ni in (111) growth sectors compared to (100). 



B. Magnetization 

Information concerning the incorporation of Ni can 
also be obtained from magnetization measurements. Be- 
tween room temperature and o:^ 100 K, the magnetic mo- 
ment of Sample B is proportional to field (Fig. 2(a)) 
with a negative slope, i.e. negative susceptibility of 
X = —2.5 X 10~^. This is close to the diamagnetic sus- 
ceptibility of diamond (x = —2.2 x 10~^) and as ex- 
pected shows no temperature dependence. Below 100 K 
an additional contribution appears, increasing in inten- 
sity as the temperature decreases. This contribution is 
attributed to a paramagnetic component. This assign- 
ment is confirmed by the good agreement (see Fig. 2(a)) 
with a Brillouin function approach. 
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FIG. 2. (Color on-line) (a) SQUID measurement of Sample 
B. The symbols show a fit using a Brillouin function as de- 
scribed in the text, (b) SQUID measurement of Sample A 
at temperatures 300, 200, 50, 10, 4, 2 K. The broken lines indi- 
cates method used to determine the spontaneous magnetiza- 
tion M{T). 



where J is the Ni quantum number, g the Lande fac- 
tor, /iB the Bohr magneton, k the Boltzmann constant, 
T the temperature, N^i the concentration of paramag- 
netic centers and B the applied magnetic field. The fit 
was performed assuming gj = 2 and leaving NNi as the 
only free parameter. A value of gJ = 2 corresponds 
to isolated nickel in the Ni2" oxidation state (3(i^), as- 
suming J ^ S ^ 1 and g = 2. A concentration of 
paramagnetic centers Nn^ = 1.5 x 10^^ cm~^ was de- 
duced. Note, a value of ^ J = 3 ( J = 3/2) or lower value 
gj = 1 ( J = 1/2) would reduce or increase N^i accord- 
ingly. Nevertheless, this concentration is in the order of 
magnitude of the expected incorporated Ni concentra- 
tion. Such a value corresponds to a Ni relative concen- 
tration of 8 X 10~^ i.e. a very diluted magnetic system 
in which weak interactions between magnetic ions and 
a paramagnetic behavior are expected. We deduce that 
nickel is incorporated as isolated, non interacting para- 
magnetic centers. 

For Sample A the magnetic moment has a very differ- 
ent behavior, as seen in Fig. 2(b). A magnetic field as 
low as 0.5 T is already sufficient to saturate the mag- 
netic moment even at 300 K, suggesting a ferromagnetic 



behavior. The linear decrease observed for higher fields 
is attributed, as before, to the diamagnetic contribution 
of the diamond matrix. Plots representing the reduced 
magnetization M(T)/Ms as a function of the reduced 
temperature T /Tc should be similar for bulk nickel and 
for this ferromagnetic sample. Under such an assump- 
tion, the temperature dependence of the experimental 
magnetization (Fig. 3) corresponds to a 500 K Curie 
temperature, 20% smaller than the one of bulk nickel 
(631 K). This suggests the presence of nickel clusters, 
with a reduced Curie temperature due to size-effects. 
The low temperature spontaneous magnetic moment of 
Ms = 7.4 X 10~^ Am^ gives a nickel concentration of 
^Ni = 3.5 X 10^^ cm~^, assuming a ^J = 2, which has 
the same order of magnitude as for Sample B. 

The strikingly different magnetic behavior of both 
samples indicates the strong influence of nitrogen on the 
incorporation of Ni in diamond. Without nitrogen {i.e. 
with a much lower N concentration than Ni concentra- 
tion), as in Sample A, Ni tends to form clusters, reducing 
the concentration of isolated nickel and the crystal is fer- 
romagnetic. With nitrogen, no ferromagnetic behavior is 
observed, the nickel atoms are diluted into the crystal as 
paramagnetic centers. This suggests a lower formation 
energy for Ni-N complexes compared to Ni defects such 
as interstitial Ni, substitutional Ni and NiV. This is in 
agreement with first principles theoretical investigation 
of nickel related complexes in diamond^ ^. 
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FIG. 3. (Color on-line) Reduced magnetic moment M{T)/Ms 
versus reduced temperature for bulk nickel (line) and Sample 
A (symbols), drawn for Tc — 631 K for nickel and Tc — 500 
K for Sample A. 



III. EFFECTIVE SPIN HAMILTONIAN 

Before presenting the magneto-photoluminescence re- 
sults we briefly outline NND effective spin Hamiltonian,^^ 
for trigonal interstitial Ni+ with the electronic configu- 
ration ZdP and effective spin S = 1/2, which is required 
to understand the data. The principle crystallographic 
orientations of the HPHT diamond crystals investigated 



here are presented schematically in Fig 4. It is possible to 
define four sets of local axes Z, X, Y (see Table I) corre- 
sponding to the four different possible (111) orientations 
of the trigonal axis of the defect. The quantization axis 
Z is parallel to the C3 axis. The perturbation of the 



magnetic field is given by the following Hamiltonian, 

AH = fiBigiBxSx + giBySy + gsBzSz), (4) 

where fiB is the Bohr magneton and the quantization 
axis is the Z trigonal axis of the center. 

The energy separation of the excited and ground state 
is large enough to ignore any interaction between them. 
The secular matrix describing the perturbation of the 
magnetic field on the ground state doublet is, 



J 



gs/^sBz + A gi/J^sBx - igi/J^sBy 

If -gsl^BBz^X qiiibBx -WiI^bBy 

2 I gil-isBx + igi/UBBv gs/^sBz - A 

gil^sBx + igiiiBBv -gs/^sBz - A 
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where A c^ 2.8 meV is the spin orbit splitting of the 
ground state and the various g terms are the effective 
Lande ^-factors. The secular matrix describing the ex- 
cited state in a magnetic field is given by, 



1 



tMb 



gIBz gfBx - igfBy 

gfBx + igfBv -glBz 



(6) 



The secular matrixes can easily be diagonalized either nu- 
merically or analytically. We have done both and verified 
that the results are identical. The analytic expressions 
for the ground and excited states are. 



Label 


Z X Y 


a 

7 

<5 


[111] [112] [110] 
[111] [1I2] [110] 
[111] [112] [110] 
[ill] [112] [110] 



TABLE L Summary of the four possible orientations of the 
local {ZXY) axis of the trigonal defect along the < 111 > 
directions. The different orientations are labeled a, /3, 7, ^. 



IV. MAGNETO-PHOTOLUMINESCENCE 
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The magnetic field lifts the degeneracy of the ground 
and excited states which split into four and two levels re- 
spectively. The energy of the possible optical transitions 
is simply the energy difference between all levels in the 
^A excited state and all levels in the ^E ground state. 
Thus, taking into account the four possible orientations 
of the defect axis (see Table I) , we expect a maximum of 
32 lines in the spectrum depending upon the orientation 
of the magnetic field. When the magnetic field is aligned 
along a symmetry axis the number of transitions with 
different energies is greatly reduced. 



Micro photoluminescence (/iPL) in dc magnetic fields 
up to 28 T and macro photoluminescence (PL) in pulsed 
magnetic field up to 56 T have been performed. For both 
experiments an optical fiber was used for the excitation 
and collection of the emission from the sample. A CW 
Ti: Sapphire laser tuned to 720 nm or a CW solid state 




FIG. 4. (Color on-line) Simplified schematic of the HTHP 
synthetic diamond crystals labeling the principal crystallo- 
graphic axes. The small (113) and (101) faces of our cubo- 
octahedral crystal have been omitted for clarity. 
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FIG. 5. (Color on-line) Typical /i-photoluminescence spectra 
of Ni incorporated in HPHT synthetic diamond for Samples A 
(at two different positions) and B at T = 4.2 K and B = T. 
All spectra are taken on the (111) face. A linear background 
has been subtracted from the rather weak PL signal of Sample 
B and the curves are offset vertically for clarity. Position 2 
is close to the bottom of the (111) face, i.e. near the seed 
crystal. See /x-PL map in Fig. 11. 




laser at 660 nm was used for the excitation. The emis- 
sion spectra have been measured using a spectrometer 
equipped with a CCD camera. For the macro PL the 
sample was placed at the end of the fiber with a diame- 
ter of 50 jam^. In jaFL measurements light was focused 
on the sample using a microscope objective coupled with 
a mono-mode fiber. The size of the laser spot on the 
sample was around Ijam^. The sample was mounted on 
piezo translation stages which allows spectrally resolved 
spatial mapping in magnetic field. The measurements 
were performed at T ~ 4 K with the magnetic field ap- 
plied parallel to (Faraday configuration) and normal to 
(Voigt configuration) the direction of propagation of the 
light. A number of different orientations of the crystal 
with respect to magnetic field have been measured on 
the (111) and (001) faces. Typical optical spectra, for 
both samples, measured at low temperature and at zero 
magnetic field are presented in Fig 5. A characteristic 
zero-phonon doublet is observed at 1.401 and 1.404 eV. 
The doublet structure originates from the A ::^ 2.8 meV 
splitting of the (^E) ground state due to a combination of 
spin orbit interaction and trigonal distortion. The asym- 
metric line shape is the result of unresolved Ni isotopic 
splitting. ^^ For sample A representative spectra recorded 
at two different locations are shown. The small shift of 
the energetic position and the slightly different doublet 
splitting observed in the spectra probably indicates indi- 
cates a different local strain. ^^ 

The intensity of the 1.4 eV emission from sample B was 
too weak to be measured in pulsed magnetic field, where 
typical integration time is of the order is 2ms. Both sam- 
ples were measured using ja-FL in dc magnetic field up to 



FIG. 6. (Color on-line) Representative T = 4.2 K magneto- 
photoluminescence spectra in the Faraday configuration of 
the 1.4 eV Ni color center in HPHT synthetic diamond, (a) 
Macro-PL with B \\ [111] collected from the (111) face, (b) 
/i-PL with B II [001] collected from the (001) face. Please 
note the different magnetic field ranges. 



28T. Selected representative macro and micro magneto- 
photoluminescence spectra for sample A are presented in 
Fig 6(a-b) respectively. The spectra were taken in the 
Faraday configuration on the (111) face with B \\ [111] 
(Fig 6(a)) and on the (001) face with B \\ [001] (Fig 6 
(b)). For both experimental configuration we observe a 
splitting on the zero phonon doublet into multiple lines 
due to the Zeeman effect. In order to compare the exper- 
imental results with the predictions of the effective spin 
Hamiltonian the energy of each transition was extracted 
from the PL spectra by fitting a Gaussian function. 



A. Macro-photoluminescence on the (111) face 

Using both Faraday and Voigt configurations as appro- 
priate, the macro-photoluminescence has been collected 
on the (111) face with the magnetic field aligned along 
the experimentally available symmetry axes of the defect 
{i.e. perpendicular to the face or parallel or perpendicu- 
lar to certain edges of the face - see Fig. 4). The energy of 
the observed optical transitions for B \\ [111], B \\ [110] 
and B \\ [112] is plotted as symbols in Fig.7(a-c). Only 



a few transitions are observed with most spectra show- 
ing only four Unes. The simphcity of the observed spectra 
immediately suggests that the defect axis is preferentially 
aligned along the [111] growth direction of the face. NND 
already reported a preferential alignment of the defect 
axis along one of the (111) directions. 

Making this assumption we fit the transitions energies 
calculated with the effective spin Hamiltonian. The con- 
figuration B II [111] for a defect aligned along [111], i.e. 
the Z axis, is a particularly easy case to fit since all the off 
diagonal terms in the secular matrices are zero and only 
the parameters gs and g^ play a role {g^ which controls 
the splitting of the excited state which is not resolved 
here is determined later from the high resolution ja-FL 
measurements). Using the NND parameters in Ref.[15] 
as a starting point we have fitted the B \\ [111] data. Sub- 
sequently, the other orientations, B \\ [110] and B \\ [112] 
were simulated to extract gi and gf. The results of such 
a fit are shown by the solid red lines in Fig.7(a-c). The 
agreement is extremely good apparently confirming the 
trigonal symmetry and preferential orientation of the de- 
fect axis along the [111] growth direction. For the [111] 
aligned defect there is no difference between the B \\ [110] 
and B \\ [112] orientations since these directions represent 
the local Y and X axes which are equivalent. This is con- 
firmed by the almost identical evolution of the observed 
transition energies for the two directions (Fig.7(b-c)). 

For B II [lIO] the green broken lines are the calculated 
transition energies for defects with their axis aligned 
along the [111] and [111] directions. Clearly these tran- 
sitions are not reproduced in the data. However, from 
this orientation we can only say that the defect is either 
aligned along the [111] or the [111] directions for which 
the projection of the magnetic field onto the local axis 
have the same magnitude so that the transitions are de- 
generate (red solid lines). It is the B \\ [112] data which 
tells us definitively that the defect is aligned along [111]; 
the projection of the magnetic field is different for all 
the other (111) directions and the calculated transitions 
for the [III] (purple dotted lines) defect orientation are 
clearly not observed. 
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FIG. 7. (Color on-line) Energy of the transitions (symbols) 
for photoluminescence collected from the (111) face of the di- 
amond crystal (Sample A) as a function of the magnetic field 
applied parallel to the [111], [lIO] and [112] directions. The 
transition energies calculated with the effective spin Hamil- 
tonian are shown as red solid lines for a defect axis aligned 
along the [111] growth direction of the face. The green broken 
lines are the calculated energies for defects aligned along the 
other (111) directions. The legend labels the Z-axis of the 
defect as in Table I. 



B. Macro-photoluminescence on the (001) face 

As before, using both Faraday and Voigt configura- 
tions, macro-photoluminescence has been collected on 
the (001) face with the magnetic field aligned along 
the experimentally available symmetry axes of the de- 
fect. The energy of the observed optical transitions for 
B II [001], B II [lIO] and 5 || [110] is plotted as sym- 
bols in Fig.8(a-c). In contrast to the (111) face, we find 
that there is no preferential orientation of the defect axis. 
The observed transitions have been simulated using the 
NND effective spin Hamiltonian. This fit was performed 
simultaneously with the fit to the data on the (111) face, 
allowing a global optimization of the parameters. The 
parameters found are summarized in Table II. 



For the B \\ [001] all the possible defect alignments are 
equivalent so that no information concerning a preferen- 
tial orientation can be obtained form this data. Never- 
theless, the predictions of the effective spin Hamiltonian 
(solid red lines) fit the data very nicely. More information 
can be gained from the B \\ [lIO] and B \\ [110] orien- 
tations. In Fig.8(b-c), the red solid lines correspond to 
defects aligned along [111] or [III] and the green dashed 
lines to defects oriented along [111] or [ill] directions. 
In the case of a preferential alignment one of the two 
B II [lIO] or B II [110] orientations would have an ex- 
tremely simple spectrum composed of only four lines. 
This is clearly not the case; all predicted transitions are 
observed for both directions. Thus, at least one of the 
[111] or [III] directions and at least one of the [111] or 



Sample A: PL collected on (001) face 
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FIG. 9. (Color on-line) Energy of the observed transitions 
(symbols) for T = 4.2 K //-photoluminescence collected on the 
(001) face with B \\ [001] on samples A and B as a function 
of the magnetic field. The red solid lines are the predictions 
of the NND effective spin Hamiltonian (all defect alignments 
are equivalent for B \\ [001]). The legend labels the Z-axis of 
the defect according to Table I. 



FIG. 8. (Color on-line) Energy of the transitions (symbols) 
for photoluminescence collected from the (001) face of the di- 
amond crystal (Sample A) as a function of the magnetic field 
applied parallel to the [001], [110] and [110] directions. The 
transition energies calculated with the effective spin Hamil- 
tonian are shown as red solid lines for a defect axis aligned 
along the [111] growth direction of the face. The green broken 
lines are the calculated energies for defects aligned along the 
other (111) directions. The legend labels the Z-axis of the 
defect as in Table I. 



[Ill] directions are occupied. The results on the (001) 
face suggest that the preferential alignment along the 
[111] direction on the (111) face may be linked to the 
growth process. 



C. /i-PL measurements 

The /i-PL technique is not as useful for determining 
the symmetry of the defect since only the Faraday con- 
figuration can be used so that the only orientation avail- 
able is with the magnetic field perpendicular to the face 
under investigation. It does however have certain advan- 
tages; (i) excitation and collection is very efficient which 
allows the investigation of samples with a low emission 



intensity and (ii) spectrally resolved maps with a spatial 
resolution c^ 1/im can be made given information con- 
cerning the homogeneity of the diamond crystal. The 
/i-PL measurements have been performed at T = 4.2 K 
using a resistive magnet in static fields up to 28 T. 

In Fig. 9(a-b) emission obtained from the (001) face of 
samples A and B is presented (symbols). The results for 
the two different samples are almost identical despite the 
approximately two orders of magnitude larger nitrogen 
concentration of sample B. As the emission intensity from 
sample B was very weak, fewer transitions are resolved 
at low magnetic field. The solid red lines are the predic- 
tions of the effective spin Hamiltonian using exactly the 
same parameters as before. As for the macro-PL mea- 
surements, the fits are very good for a data set of much 
higher quality further confirming the trigonal symmetry 
of the defect. For the (001) face with B \\ [001] all the 
defect orientations are equivalent so that no information 
concerning defect alignment can be extracted. The nickel 
concentration in the two samples is approximately equal 
so that the results suggest that while excess nitrogen re- 
duces considerably the intensity of the PL emission, pre- 
sumably by forming other complexes with Ni, some of 
the Ni nevertheless enters the host lattice as 1.4 eV Ni 
color centers. 

Fig. lO(a-b) shows /i-PL measurements obtained at two 
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FIG. 11. (Color on-line) Spatial PL map of the (111) face 
measured at T = 4.2 K and B = T showing the integrated 
intensity of the high energy component of the 1.4 eV doublet. 
The dashed line indicates the bottom of the (111) face. Inten- 
sity decreases rapidly below this line due to the loss of focus. 
The positions at which the /i-PL spectra in Fig. 5 were taken 
are indicated. 



FIG. 10. (Color on-line) (a-b) Energy of the observed transi- 
tions (symbols) for T = 4.2 K /x-photoluminescence collected 
at two different positions on the (111) face with B \\ [111] 
on sample A as a function of the magnetic field. The pre- 
dictions of the NND effective spin Hamiltonian are shown 
for a defect preferentially aligned along [111] (red solid lines) 
and the other three equivalent (111) directions (green broken 
lines). Note that at Position 2 the defects are not preferen- 
tially aligned along the growth direction and fitting requires 
a radically different set of parameters, (c) Splitting of the ex- 
cited state; the high (circles) and low (triangles) energy tran- 
sitions in (a) , from which the average energy of the transition 
has been subtracted, are plotted versus magnetic field. The 
solid lines are a least squares fit to extract ^| = 0.16 it 0.03. 



different positions on the (111) face of sample A with 
B II [111]. We have performed a full map of the face and 
position 1 in Fig. 10(a) is representative of most of the 
(111) face. The behavior is identical to that observed in 
macro-PL. The predictions of the effective spin Hamil- 
tonian for a defect preferentially aligned along the [111] 
direction are shown by the red solid lines and the other 
orientations by the green broken lines. The data is well 
fitted by the preferentially aligned defect scenario. As 
emission from position 1 is representative of the (111) 
face this allows us to conclude that the defect is prefer- 
entially aligned along [111] growth direction over most of 



the (111) face. The splitting of the excited state is also 
clearly seen in this high resolution low noise data taken 
in static magnetic fields. See e.g. the splitting of the 
lowest and highest energy transition in Fig. 10(a) above 
15 T. Fitting the Hamiltonian to this splitting it is pos- 
sible to extract a refined value of g^ which describes the 
small spin splitting of the excited state with the field ap- 
plied along the Z quantization axis. In Fig. 10(c) the 
excited state splitting is shown; we plot the energy of the 
split transitions after subtracting the average energy to 
remove the background. A least squares fit to both data 
sets (solid lines) gives g^ = 0.16 ±0.03. The complete set 
of parameters for the effective spin Hamiltonian is given 
in Table 11. The values proposed by NND are shown for 
comparison. The differences are small showing both the 
reproducibility of the results between different samples 
and the remarkable job performed by NND from limited 
PL data in relatively low magnetic fields B < 6 T. 

It is however possible, towards the bottom of the (111) 
face i.e when going far away from the [lIO] top edge 
(see Fig. 4), to find very different spectra as shown in 
Fig. 10(b). The spectra measured at position 2 contain 
many more lines indicating that the defect is not prefer- 
entially aligned along the growth direction. Such spectra 
occur in areas of the (111) with increased intensity of 
emission in zero magnetic field (see Fig. 11). In addi- 
tion a radically different parameter set for the effective 
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Standard parameters (this work) 


2.3 ±0.05 


2.0 ±0.05 


1.7 ±0.05 


0.16 ±0.03 2.4 ±0.05 


2.800 


Defect misaligned on (111) 


face 


3.5 ±0.1 


3.0 ±0.1 


1.7 ±0.1 


(0.16) 


2.3 ±0.1 


2.874 



TABLE IL Parameters of the effective ^S = 1/2 spin Hamiltonian found here and in previous work. The standard parameters 
fit the majority of the present data. The misahgned parameter set corresponds to the special case of defects on the (111) face 
which do not have their quantization axis preferentially aligned along the [111] growth directions. The larger estimated error 
for the misaligned defects is due to the restricted data set which is limited to the B \\ [111] configuration on the (111) face. 
Note, the NND values are twice those given in Ref.[15] due to their implicit use of an 5* = 1 Hamiltonian when writing the 
secular equation. 



spin Hamiltonian is required to fit the data. We have 
fitted the data assuming that the defect orientation is 
not along the [111] growth direction (green broken lines). 
For B II [111] the three "misaligned" defect orientations 
are all equivalent. The extracted parameters are given in 
Table H. The fit is almost perfect confirming the trigo- 
nal symmetry. The red solid lines are the prediction for 
a defect aligned along the [111] growth direction. From 
the data above 20 T it appears that these transitions 
are completely absent from the spectra. This suggests 
that the defect may even be preferentially "not aligned" 
with the growth direction in this region of the face which 
is close to the seed crystal. Moreover, the parameters 
required to fit are considerably different from those ob- 
tained previously. In particular, the values of gs and 
g^ are ~ 50% larger reflecting the much larger Zeeman 
splitting of the ground state even though the magnetic 
field is not aligned along the Z quantization axis. This 
suggests that the local environment (trigonal distortion) 
of the Ni center is significantly different in regions of the 
(111) face where it does not preferentially align along the 
[111] growth direction. 



V. DISCUSSION 

The 1.4 eV Ni center investigated here has been unam- 
biguously identified with the NIRIM-2 ESR line in HPHT 
diamond. ^^'^^ NIRIM-2 has a large angular dependence 
of the magnetic field position of the ESR lines which is 
consistent with the center having trigonal symmetry. ^^ 
ESR gives very precise values for the g-factors but probes 
only the splitting of the ground state. While magnetic 
circular dichroism (MCDA) is a less precise technique, it 
is extremely useful since it gives access to the ^-factors 
of the excited state. ^^ The measured anisotropy of the g- 
factor can provide a crucial test for a given microscopic 
model. In Fig. 12 we show the ground state splitting (ef- 
fective ^-factor) obtained from the ESR data of Isoya et 
al}^ together with the MCDA results of Mason et al}^ 
for the excited state splitting. The expected angular de- 



pendence of the ground state splitting using the effective 
spin Hamiltonian is also plotted for the ground (solid 
lines) and excited state (broken line) as a function of an- 
gle when the magnetic field is rotated around the [110] 
direction, e.g. from B \\ [001] to B \\ [112]. The pre- 
dictions of the effective spin Hamiltonian are in excellent 
agreement with the ESR and MCDA results for both the 
ground and excited states. 

In order to compare the ^-factors of the effective spin 
Hamiltonian with the ESR and MCDA results, it is con- 
venient to define effective ^-factors, g\\ and g^ corre- 
sponding to a magnetic field aligned parallel and per- 
pendicular to the Z quantization axis of the defect. For 
the ground state we have a large g\\ = gs and g± = 0. 
The magnetic field component which is not along Z does 
not generate a Zeeman splitting of the ground state. The 
Bx and By components only change the splitting of the 
zero field doublet as can be seen from Eq. 7. The situ- 
ation for the excited state is very different with a small 
^11 = ^3 and a large g± = gf. As pointed out by Mason et 
al.^ the very small value of ^y ~ is a most unexpected 
result for an orbital singlet state. ^^ The high magnetic 
field PL results presented here have allowed a rather pre- 
cise determination of the ^-factors which motivates us to 
revisit Mason et al.^s quantitative comparison with the g- 
factors derived from the crystal field theory of interstitial 
3d^ Ni+. 



A. Crystal field model 

It has been proposed that the Hamiltonian used for 
neutral substitutional vanadium in SiC^^ can also be ap- 
plied to the case of interstitial Ni+,^^ 



H = Ht^ + Hc:^^ + Hs.o + Hz 



(8) 



where Ht^ contains the cubic crystal field term, i^Cs^ the 
trigonal crystal field, Hs.o the spin-orbit interaction and 
finally Hz includes the Zeeman interaction. The result- 
ing level structure is shown in Fig. 13. The parameters 
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FIG. 12. (Color on-line) The circles show the angular depen- 
dence of the line position (effective ^-factor) of the NIRIM-2 
center, taken from Ref.[13], with the magnetic field rotated 
about the [110] axis at T = 4 K and 5 :^ 0.5 T. The tri- 
angles are the splitting of the excited state measured using 
MCDA taken from Ref.[17]. Also shown is the predicted an- 
gular dependence of the ground state splitting (solid lines) 
and excited state splitting (broken lines) , calculated using the 
effective spin Hamiltonian with the parameters found in this 
work, and given in Table II. The curves are labeled according 
to the alignment of the defect Z quantization axis. 



Ac , iC, C^ cind k are the cubic crystal- field splitting, one- 
third of the trigonal crystal-field splitting within the ^T2 
state, the spin-orbit parameter, and an orbital reduction 
factor, respectively {K' ^ ("', and k' are the corresponding 
quantities for the ^E state). The sign of the ^T2 trigo- 
nal crystal field parameter, iC < 0, is chosen so that the 
Ai level of the excited state is lowest in agreement with 
the uniaxial stress measurements^^ and the spin-orbit pa- 
rameters ( and C' are both negative for the single (P hole. 
The ^E ground state is split into two Kramers doublets 
by a combination of spin-orbit coupling and trigonal dis- 
tortion. 



A = ^(r4)-^(r5,6) 



4CK' 



(9) 



where the cubic crystal field term Ac — 1.4 eV. From the 
measured splitting A = 2.8 meV we obtain the product 
CK' = 0.98 X 10-^ eV^ 

The ^-factors of the ground state doublet are for the 
Fs 6 state, -^'^ 



9\\= 93 = ^ 



4kX' ^k'K' 



FIG. 13. (Color on-hne) Energy levels of the '^D state of 3d^ 
interstitial Ni^ showing the subsequent effects of the cubic 
{Td) and trigonal (Cs^) crystals fields and spin-orbit inter- 
action.^^ The broken arrows indicate the observed zero field 
optical transitions. 



and for the F4 state 

^11 = ^3 = 2 

4.k'C' 
91- = -1^' 



4kX' ^k'K' 



(11) 



(10) 



As Ac ^ K' ^ C,' the F4 state has ^x ^ so that the ef- 
fective spin Hamiltonian (in which de facto g±_ = 0) is a 
reasonable approximation. Using the values of ^3 = 2.3 
and ^3 = 2.0 found in this work, solving Eqs.( 10-11) im- 
mediately gives Ak'C l^c = Sk'K^Ac = -0.15 so that 
g± = —0.15. Although small, the splitting of the up- 
per doublet state should be experimentally observable at 
magnetic fields above 15 T with B applied perpendicular 
to the Z defect axis on the [111] face. The high resolu- 
tion /i-PL measurements, with its narrower line widths 
and better signal/noise in static magnetic fields, clearly 
resolved the similarly small splitting of the excited state 
with B \\ Z. Unfortunately, only the Faraday configura- 
tion is possible when using the /i-PL system so that such 
a measurement is excluded for the present. Making the 
approximation that the orbital reduction factor /c' ~ 1 
gives K' = C/2 = 0.0262 eV. As pointed out by Mason 
et al.^ while such values of (' and K' are reasonable, the 
product ('K' = 1.4 X 10~^ eV^ is too large to be compat- 
ible with the CK' - 1 x 10"^ eV^ obtained from Eq.(9) 
knowing the spin-orbit splitting A = 2.8 meV of the zero 
field doublet. Using an orbital reduction factor k' < 1 
only makes the situation worse. 

The ^-factors of the excited state (lower F4 of ^T2) are 
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given by. 
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VI. CONCLUSION 



g±= gl = 2a^ - 2V2kab, 



(12) 



where, 



a = cos(7), 
b= -sin(7), 

C/V2 



3i4:/2-C/4' 



(13) 



The cos and sin terms give an additional constraint 
that a^ + 6^ = 1, so that the orbital reduction factor k 
cannot be freely chosen. Using the values of ^1 = 0.16 
and gf = 2.4 found in this work, a solution to Eq.(12) 
exists with k = 0.79, a = 0.82 and b = -0.574. From 
Eq.(13) this implies that K :^ C/3.5. Neglecting cova- 
lency effects, which are expected to be small in diamond, 
X^ y giving K ^ 0.015 eV. 

Thus the only apparent short coming of the crystal 
field model is its inability to correctly predict the zero 
magnetic field doublet splitting. It has been suggested 
that this might be explained by corrections to the split- 
ting of the '^E state which can arise due to a mixing with 
the 4p states of the Ni+.^^ Such a mechanism has also 
been proposed for 3dP Cu^+ in II- VI compounds. ^^'^^ 

A crucial test of the crystal field model is provided 
by defects with a quantization axis which is not pref- 
erentially aligned: In the /i-PL data it is possible to 
find regions of the (111) face, towards the bottom of 
the sample, i.e. near the seed crystal, where the de- 
fects are not preferentially aligned along the [111] growth 
direction. Such defects have an almost identical emis- 
sion energy (see Fig. 5) and an almost identical splitting 
(A = 2.874 meV) of the zero magnetic field doublet. Sur- 
prisingly, the ^-factors for the trigonal axis of such defects 
are markedly different (see Table. II). Using the values 
found in experiment, gs = 3.5 and g'^ = 3.0 in Eqs.(10- 
11), gives 4A:'C7Ac = -1.25 and Sk'K^Ac = -0.25 so 
that K^ = CV^O- Assuming k^ ^ 0^ the spin-orbit cou- 
pling term (' c^ 0.44 eV and trigonal crystal field term 
K' :^ 0.044 eV. The value of (' is an unreasonable 6 times 
larger that the value of 0.075 eV for a free Ni+ ion.^^'^^ 
Moreover, the product CK' !:^ 20 x 10~^ eV^ overesti- 
mates the zero field doublet splitting (Eq.(9)) by over an 
order of magnitude. At the same time experiment indi- 
cates that the splitting of the zero field doublet is almost 
unchanged. This forces us to conclude that the crystal 
field model does not provide a correct quantitative de- 
scription of the 1.4 eV Ni color center in diamond. 



We have investigated two different boron-free HTHP 
synthetic diamond crystals. Our high magnetic field 
photoluminescence results are perfectly described by the 
NND effective spin Hamiltonian, unequivocally demon- 
strating the trigonal symmetry of the 1.4 eV Ni color cen- 
ter in diamond. Both samples investigated have a similar 
Ni content (c:^ 1 x 10^^ cm~^) but radically different con- 
centrations of N. Sample B has a similar concentrations 
of Ni and N, while sample A has roughly two orders of 
magnitude less nitrogen. Magnetization measurements 
show that Ni is predominantly incorporated into Sample 
B as isolated Ni complexes and into Sample A as nm 
size Ni clusters. Nevertheless, both samples exhibit the 
characteristic 1.4 eV emission doublet associated with an 
isolated Ni complex. Despite the significantly lower con- 
centration of isolated Ni, Sample A shows much stronger 
1.4 eV emission. This suggests that the presence of N 
does not necessarily favor the formation of the 1.4 eV 
Ni color center; the concentration of isolated Ni-N com- 
plexes is almost certainly higher in Sample B, while the 
1.4 eV PL signal is much weaker, suggesting that the 
presence of N may actually impede the formation of this 
particular Ni color center. As nitrogen is usually a donor 
in diamond, a possible mechanism could be the transfer 
of a donor electron to the Ni+ ion reducing the number 
of optically active centers. ^^'^^ 

The magneto-PL presented here, together with pre- 
viously published magneto-PL and PL under uniaxial 
stress, ^^'^^ ESR^^ and MCDA^'^ are ah consistent with 
an interstitial 3cfi Ni+ with spin 6* = 1/2 and a large trig- 
onal distortion {Csy symmetry) due to a displacement of 
the Ni along a < 111 > direction. The exact nature of the 
complex nevertheless remains unknown. A complex in- 
volving an additional, transition metal impurity, dopants 
such as nitrogen or boron, or a vacancy or divacancy, is 
required to produce the trigonal distortion. On the basis 
of the ESR results Isoya et al. suggested an intersti- 
tial A^i+-vacancy complex. However, this configuration 
has been shown to be unstable^^ with the interstitial Ni 
moving towards vacancy i.e. to the substitutional site. 
First principles calculations suggests that that a com- 
plex involving boron and substitutional nickel fulfills all 
the necessary requirements.^^ However, this assignment 
seems to be unlikely here as there was no boron present 
in the melt during the growth of our samples. Moreover, 
recent density functional theory (DFT) calculations cou- 
pled with X-ray absorption spectroscopy measurements 
performed on Sample B, suggest the incorporation of 
nickel as a divacancy complex, in which interstitial Ni 
is placed at the midpoint between two vacancies. ^^ How- 
ever, the first principles calculations of Larico et al.^^ 
suggest that the (VNiV)+ complex has C2/1 symmetry 
with the Ni related electronic states resonant and inert 
inside the valence band. First principles calculations of 
defects in diamond are in general complicated due to the 
need to take into account possible relaxations of neigh- 
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boring carbon atoms which can alter both the energy 
and the symmetry of the center. Under such conditions 
the microscopic model for the 1.4 eV Ni color center in 
diamond should still be considered as an open question. 
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